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1. Introduction 

In high energy hadron-hadron collisions, jet production is the dominant high transverse mo- 
mentum process. It results from a scattering of a coloured parton off another coloured parton, a 
process whose cross section is described in the following factorised form, 

d<T = E / ^^fi{^^F)fj{^2^F)^dij{a s {lXR),lXR,lX F ). (1.1) 
ijJ Si S2 

In ( |1.1D the probability of finding a parton of type i in the proton carrying a momentum fraction 
is described by the parton distribution function fi(^ ,/j.p)dt l and the parton-level scattering cross 
section da, for parton i to scatter off parton j normalised to the hadron-hadron flux 1 is summed 
over the possible parton types i and j. As usual }Xr and jj,p are the renormalisation and factorisation 
scales. 

The inclusive jet pp — > j + X and dijet pp —> jj + X cross sections have been studied at the 
TEVATRON ^ and at the LHC [Q] where the advantage of very high statistics and steady im- 
proved control over systematic uncertainties ^ allows these measurements to be performed with 
a few percent accuracy. This high quality jet data which probes the dynamic of the hard scattering 
simultaneously at a wide range of momentum transfers thus contains precision information on the 
structure of the proton and on the QCD coupling constant a s . However, extraction of these funda- 
mental parameters from the data depends on reliable theoretical calculations which should match 
the precision of the experimental measurement. 

The uncertainties related to the theoretical prediction may be improved by including higher or- 
der perturbative corrections which have the effect of (a) reducing the renormalisation/factorisation 
scale dependence of the cross section and (b) improving the matching of the parton level event 
topology with the experimentally observed hadronic final state [01]. In this way the partonic cross 
section da, in ( fLl| ) has the perturbative expansion 

da,- = doff + (H) daf o + 2 dd NN w + ^ (a 3 } (L2) 

where the next-to-leading order (NLO) and next-to-next-to-leading order (NNLO) strong correc- 
tions are identified. The NLO corrections for single inclusive jet and dijet distributions have been 
computed by several groups []|] and have been successfully compared with data from the TEVA- 
TRON [jl]| and LHC In this talk we are mainly concerned with the pure gluon channel contri- 
bution to the NNLO corrections to dijet production where the resulting theoretical uncertainty is 
estimated to be at the few per-cent level [Q]. 



2. NNLO calculations 

At NNLO, there are three distinct contributions due to double-real radiation da§^ LO , mixed 
real- virtual radiation da^J /L0 and double- virtual radiation da^ LO , that are given by 

dO NNL0 = d&NNW+ f d &MLO+ f dd NNLO C 2 " 1 ) 

■Jd<S> m+2 Jd®,„ +1 Jd® m 

'The partonic cross section normalised to the parton-parton flux is obtained by absorbing the inverse factors of |] 
and £2 into drj/y. 
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where the integration is over the appropriate A^-particle final state subject to the constraint that 
precisely ra-jets are observed, 

/ = / dfctf/W. (2.2) 
For dijet production in the pure gluon channel these contributions read, 

d ^NNLO = N J d ^4iP3,P4,P5,P6;PUP2)\^ggKgggg\ 2 J2^(j>3,P4,P5,P6) 
d &MLO = N J d ^3(P3,/H,/>5;/?l,/>2)(^|^ 

where denotes the /-loop amplitude and the jet function J$ defines the procedure for building 
ra-jets from «-final state partons. 

Individually the partem channels defined above are all separately infrared divergent although, 
after renormalisation and factorisation of initial state singularities, their sum is finite [|, 0. The sin- 
gularities in the double-real channel arise in the single and double unresolved regions of the phase 
space which physically correspond to two gluons becoming simultaneously soft and/or collinear. 
In these regions of the phase space the tree-level matrix element diverges in universal and process 
independent ways The real- virtual contribution contains, on the one hand, explicit infrared di- 
vergences coming from integrating over the loop momentum and, on the other hand, implicit poles 
in the regions of the phase space where one of the final state partons becomes unresolved. This cor- 
responds to the soft and collinear limits of the one-loop amplitude which again have universal and 
process independent limits 0. Finally the double- virtual contribution contains the renormalised 
two-loop corrections containing explicit infrared poles coming from integrating over the loop mo- 



menta [10] 



We note that explicit expressions for the interference of the four-gluon tree-level and two- 
loop amplitudes are available in Refs. JTT[], while the self interference of the four-gluon one-loop 



amplitude is given in [12]. The one-loop helicity amplitudes for the five gluon amplitude are given 



in Q13| ] and the double-real six-gluon matrix elements were derived in Q14|]. All matrix elements 
relevant for the computation of gg — >• gg at NNLO are therefore known. 

The bottleneck for the computation of jet observables at NNLO in this case is that the real 
and virtual contributions have a different number of final state partons and require the phase space 
integrations to be done separately. As a result one is left with the task of finding a systematic 
procedure to extract the infrared soft and collinear singularities that are present in the double-real 
radiation and mixed real-virtual contributions at NNLO before a finite physical prediction can be 
obtained. 

Subtraction schemes are a well established solution to this problem and several methods 
for systematically constructing general subtraction terms have been proposed in the literature at 



NLO [pp and, with various degrees of sophistication, at NNLO [17, 23-E5h. Another NNLO 



subtraction scheme has been proposed in [18]. It is not a general subtraction scheme, but can 
nevertheless deal with an entire class of processes, those without coloured final states, in hadron- 
hadron collisions. 
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In addition, there are purely numerical methods [19-21] such as the completely independent 
approach called sector decomposition, which avoids the need for analytical integrations and which 
has been developed for virtual fli9| ] and real radiation [ |20| ] corrections to NNLO. Finally, a new 
numerical method combining the ideas of subtraction and sector decomposition has been proposed 



and developed in [22]. We will follow the NNLO antenna subtraction method which was proposed 



in [23]. 



3. Antenna subtraction 

The subtraction approach is based on the universal factorisation of the QCD matrix elements 
(as described in the previous section) and of the phase space in the singular regions. 

The factorisation property of the matrix elements allows us to find simple counterterms for the 
real-radiation channels which reproduce the limits of the amplitudes in the different singular kine- 
matic configurations and therefore regulate the numerical integration of the real contribution. The 
phase space factorisation property on the other hand allows us to carry out the analytic integration 
of the countertem such that it can be added to the virtual contribution in integrated form leading to 
an analytic pole cancellation between the two. This analytic integration is fully inclusive over the 
unresolved configurations and can be carried out as it is not subject to any jet defining cuts, and 
does not depend on the details of the process under consideration. The antenna subtraction scheme 
is then by construction a general subtraction scheme and, as we will describe below, can be applied 
to any NNLO QCD calculations for exclusive jet observables at hadron colliders. 

Suppressing the labels of the partons in the initial state of the hard scattering and introducing 



subtraction terms, equation (gjj) can be rewritten as [ |23| ] : 



d&NNT.n = I (dOMMm-dduurn) + / da, 5 



NNLO - / {UO NNL0 -UO NNL0 ) i- / UOJVJVLO 

JdQ> m+ 2 Jd<&„, + 2 

+ / (d< t0 - da v NNLO ) + / da v N s NL0 + / daZw 

d&NNLO+ / dOxNLO- C 3 - 1 ) 

Here, ddf lNL0 denotes the subtraction term for the (;n + 2)-parton final state which behaves like the 
double-real radiation contribution da^ LO in all singular limits. Likewise, da^ L0 is the one-loop 
virtual subtraction term coinciding with the one-loop (m+ 1) -final state da^ LO in all singular 
limits. The two-loop correction to the m-parton final state is denoted by dd^ L0 . In addition, as 
there are partons in the initial state, there are also two mass factorisation contributions, da^ L l and 
^nnlo' f° r tne ( m + 1)" an d m-particle final states respectively. 

In order to construct a numerical implementation of the NNLO cross section, the various 
contributions must be reorganised according to the number of final state particles, 



dGNNLO — / [d&NNLO ~ ^Nl> 

Jfr& m+1 

+ / [dOmw-doL 

Jd<S> m+l 

+ I Wnlo-Mnnlo], (3-2) 

Jd<I>,„ 



>NNLO\ 
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where the terms in each of the square brackets are finite and well behaved in the infrared singular 
regions. More precisely, 

ao NNLO — ao NNLO~ J^ ao NNLO~ a °NNLOi ^.J) 

&°NNLO = ~ J d °NNLO ~ J d °NNLO ~ dd NNLO- ( 3 " 4 ) 

where the double-real integrated subtraction term gives contributions to both (m + l)- and m-parton 
final states, 

f da s mw =( Ua^ w + f Ud s / NW . (3.5) 
The derivation and specific forms of the dd^ NL0 and dd^ L0 counterterms for the case of 



coloured particles in the final state has been discussed in [23]. The extension to the case of coloured 



particles in the initial state relevant for studying processes at hadron colliders is described in [24, 



25] 



In a brief description, the subtraction terms are constructed from products of antenna functions 
with reduced matrix elements (with fewer final state partons than the original matrix element). 
The relevant antennae are determined by both the external state and the pair of hard partons it 
collapses to. In general we denote the antenna function as X. For antennae that collapse onto a 
hard quark-antiquark pair, X = A for qgq. Similarly, for quark-gluon antenna, we have X = D for 
qgg and X = E for qqq final states. Finally, we characterise the gluon-gluon antennae as X = F 
for ggg, X = G for gqq final states. In all cases the antenna functions are derived from physical 



matrix elements: the quark-antiquark antenna functions from y* — > qq + (partons) [26], the quark 



gluon antenna functions from % — > g + (partons) [ 27 ] and the gluon-gluon antenna functions from 
H — > (partons) [|28|]. The full subtraction term is obtained by summing over all antennae required 
for the problem under consideration. In the most general case (two partons in the initial state, and 
two or more hard partons in the final state), this sum includes final-final, initial-final and initial- 
initial antennae. The antenna subtraction formalism has also been extended to deal with hadronic 
observables involving massive fermions at NNLO [ |2^ ] . 

A key element in any subtraction method, is the ability to add back the subtraction term in- 
tegrated over the unresolved phase space. The subtraction terms are then integrated over a phase 
space which is factorised into an antenna phase space (involving all unresolved partons and the two 
radiators) multiplied with a reduced phase space (where the momenta of radiators and unesolved 
radiation are replaced by two redefined momenta). In the antenna approach, this integration needs 
to be performed once for each antenna and the analytic integration of those is performed over the 
antenna phase space only. All the necessary integrals relevant for jet cross sections at NNLO at 
hadron-hadron colliders have been computed and are summarised in Table 1 . 



4. Dijet production at NNLO at hadron colliders 

In this section we will discuss the ongoing calculation of dijet production at NNLO at hadron 
colliders using the antenna subtraction scheme. 
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Table 1: Integrated massless antenna functions 
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4.1 Double-real contribution 



In a previous paper [1240, the subtraction term ddtj, NLO in (|3.2|), corresponding to the leading 



colour pure gluon contribution to dijet production at hadron colliders was derived. The subtraction 
term was shown to reproduce the singular behaviour present in da^ L0 in all of the single and 
double unresolved limits. One example of this can be seen in fig. [j](b) where we generated 10000 
random double soft phase space points and show the distribution of the ratio between the matrix 
element and the subtraction term Rrr. The three colours represent different values of x = (s — Sij)/s 
[x = 10~ 4 (red), x = 10~ 5 (green), x = 10~ 6 (blue)] and we can see that for smaller values of x we 
go closer to the singular region and the distribution peaks more sharply around unity. 



double soft limit for gg .gggg 

#PS points=10000 
1487 outside the plot 

317 outside the plot 
59 outside the plot 



x=1(H ™ 

x=10 "fi ™ 
x=10" 6 H 

X=(S-Sy)/S 



0.99997 0.99998 0.99999 



1.00001 1.00002 1.00003 



(a) 



(b) 



Figure 1: (a) Example configuration of a double soft event in the double-real channel with Sjj 
Distribution of Rrr for 10000 double soft phase space points. 



i s n = s. (b) 



We also analysed the regions of phase space corresponding to the final and initial state single 
collinear singularities. There we saw that the subtraction term, which is based on azimuthally 
averaged antenna functions does not accurately describe the azimuthal correlations present in the 
matrix elements and antenna functions when a gluon splits into two collinear gluons. Nevertheless, 
the azimuthal terms coming from the single collinear limits do vanish after an azimuthal integration 
over the unresolved phase space with their angular dependence being of the form Acos(20 + a), 
where (j> is the azimuthal angle of the collinear system. Therefore by combining two phase space 
points with azimuthal angles <p and + § [ ^4] , 35 ] and all other coordinates equal, the azimuthal 
correlations drop out since the correlations present in the matrix element at the rotated point cancel 
precisely the azimuthal correlations of the un-rotated point. 
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#phase space points = 10000 
5 outside the plot 
outside the plot 
outside the plot 



10000 | 1 1 

flphase space points = 10000 
221 outside the plot 
28 outside the plot 
8000 - 1 outside the plot 



1 1.001 1.002 1.003 1.004 1.005 



(a) 



(b) 



Figure 2: Distribution of Rrv for 10000 single collinear phase space point pairs where the pair of phase 
space points is related by an azimuthal rotation ofn/2 about the collinear direction, (a) Final state collinear 
singularity (b) Initial state collinear singularity 



4.2 Real-virtual contribution 



In a subsequent paper [25] we derived the subtraction term da^ NLO in (3^2) corresponding 



to the leading colour pure gluon contribution to dijet production at hadron colliders. In [25] we 
established which contributions from the double-real subtraction term have to be added back in 
integrated form at the (m+ 1) and ra-parton level. By construction the counterterm removes the 
explicit infrared poles present on the one-loop amplitude, as well as the implicit singularities that 
occur in the soft and collinear limits. The £-poles present in the real-virtual contribution are an- 
alytically cancelled by the e-poles in the subtraction term rendering the real-virtual contribution 
locally finite over the whole of phase space. 

We also tested how well the real-virtual subtraction term dd^ NL0 approaches the real-virtual 
contribution dd§^ jL0 in all single unresolved regions of the phase space so that their difference can 
be integrated numerically in four dimensions. This is shown in fig. || for the single collinear case. 
We noted [|5|] that in the real-virtual channel the magnitude of the angular correlations is much 
smaller than in the double-real case. Nevertheless the combination of two phase space points with 
azimuthal angles and + f can be used to provide an improvement in the convergence of the 
subtraction term. 

4.3 Double-virtual contribution 



For the double-virtual contribution we are concerned with the NNLO contribution coming 
from two-loop processes, dd^ L0 , the remaining integrated subtraction terms and NNLO mass 



factorisation contributions that are collectively denoted da^ NLO in (3^4). We note that this piece 



includes the full integrated real- virtual subtraction term do^ L0 derived in [25] as well as integrated 



contributions from the double-real subtraction term derived in [24]. This involves the integrated 



initial-initial four parton antenna functions which were recently computed in [p3|], the integrated 

I We 



one-loop three parton antennae [ 34 ] and convolutions of integrated three parton antennae 
have verified that explicit poles from the double- virtual matrix elements cancel analytically against 
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poles from this subtraction term. More explicitly [p 



@>oles \ dd^ L0 



NNLO 



cMF. 



2 \ 



NNLO 



(4.1) 



where the operator J^oles selects the singular contributions containing £ -poles that start at ff(e 



-4\ 



This gives us strong confidence on the derivation of the subtraction terms in [|24[ [25J] and on the 



correctness of the results presented in []33[]. In addition, it tells us that the infrared singular structure 
of the pure gluon scattering matrix elements at NNLO is captured analytically by the antenna 
subtraction method in a systematic and accurate manner. 



5. Conclusions 



In this talk we discussed the generalisation of the antenna subtraction method for the calcu- 
lation of NNLO QCD corrections for exclusive collider observables with coloured partons in the 
initial state. We considered the pure gluon channel contributing to the dijet cross section at NNLO 
and showed that the explicit £ -poles in real-virtual and double-virtual contributions are analytically 
cancelled by the £-poles in the subtraction terms. The finite remainders in the double-real, real- 
virtual and virtual-virtual channels can be evaluated numerically in four dimensions by generating 
appropriate weighted partonic events. 

Future steps include the construction of a numerical program to compute NNLO QCD correc- 
tions to dijet production in hadron-hadron collisions with the inclusion of the remaining channels 
involving quark initiated processes [36]. 
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